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Let n„ be the set of convex polygonal lines F with vertices on 
1.^ , Z\_ and fixed endpoints — (0, 0) and n = (ni, 712). We are concerned 

j-^ • with the limit shape, as n — ^ 00, of "typical" F £ i7„. with respect 

to a parametric family of probability measures {P^, < r < 00} on 
fT^ ■ TTn, including the uniform distribution (r = 1) for which the limit 

shape was found in the early 1990s independently by A. M. Vershik, 
I. Barany and Ya. G. Sinai. We show that, in fact, the limit shape 
is universal in the class {Pn}, even though P^ {r ^ 1) and P^ 
are asymptotically singular. Measures P^ are constructed, following 
ri . Sinai's approach, as conditional distributions Qz(-|-ffn), where Ql 

are suitable product measures on the space 77 = Un77„, depending 
on an auxiliary "free" parameter z = (21,22). The transition from 
{n,Q'^) to {FlnjPn) is based on the asymptotics of the probability 
Ql{n„), furnished by a certain two-dimensional local limit theorem. 
*vj , The proofs involve subtle analytical tools including the Mobius inver- 

^ ■ sion formula and properties of zeroes of the Riemann zeta function. 

(N 

00 

^ ' 1. Introduction. 

rn 

r^ I 1.1. Background: the limit shape. In this paper, a convex lattice polyg- 

QQ ' onal line i^ is a piecewise linear path on the plane, starting at the origin 

O ■ = (0,0), with vertices on the integer lattice Z^ := {{i,j) S Z^ : «,J > 0}, 

and such that the inclination of its consecutive edges strictly increases stay- 
ing between and 7r/2. Let 77 be the set of all convex lattice polygonal lines 
with finitely many edges, and denote by 7T„ C 77 the subset of polygonal 
5^ i lines 7~' G 77 whose right endpoint C = Cr is fixed at n = (ni, 712) S Z^ . 

We are concerned with the problem of limit shape of "typical" 7^ G 77.„ , as 
n — )• 00, with respect to some probability measure P„ on 77„. Here the "limit 
shape" is understood as a planar curve 7* such that, with overwhelming P„- 
probability for large enough n, properly scaled polygonal lines Fn = Sn{r) 
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lie within an arbitrarily small neighborhood of 7*. More precisely, for any 
e > it should hold that 

(1.1) lim P„{d(f„,7*)<e} = l, 

n— >oo 

where d{-,-) is some metric on the path space — for instance, induced by the 
Hausdorff distance between compact sets (in M?), 



(1.2) d-}{(^, B) := max < maxmin Ix — yl, maxmin \x — y\ 

where | • | is the Euclidean vector norm. 

Of course, the limit shape and its very existence may depend on the 
probability law P„. With respect to the uniform distribution on 77„, the 
problem was solved independently by Vershik [32], Barany [3] and Sinai 
[29], who showed that, under the scaling 5'.„ : (xi,X2) 1— )• (xi/ni, X2/n2), the 
limit shape 7* is given by a parabola arc defined by the Cartesian equation 

(1.3) Vl - 3;i + Vx2 = 1, 0<xi,X2<l. 

More precisely (cf. (1.1)), if n = (71-1,722) — )■ 00 so that n2/ni — )• c G (0,oo) 
then, for any e > 0, 

(1.4) um *ir<=.n„:Mr„.,')<^) ^ , 

' »^~ #(ir„) 

(Here and in what follows, ij^{-) denotes the number of elements in a set.) 

The proofs in papers [32, 3] involved a blend of combinatorial, variational 
and geometric arguments and were based on a direct analysis of the corre- 
sponding generating function via a multivariate saddle-point method for a 
Cauchy integral [32] or a suitable Tauberian theorem [3]. Extending some of 
these ideas and using large deviations techniques, Vershik and Zeitouni [37] 
developed a systematic approach to the limit shape problem for the uniform 
measure on more general ensembles of convex lattice polygonal lines with 
various geometric restrictions. 

Sinai [29] proposed an alternative, probabilistic method essentially based 
on randomization of the right endpoint of the polygonal line F G 7T„ ; we will 
comment more on this approach in Section 1.3. Let us point out that the 
paper [29] contained the basic ideas but only sketches of the proofs. Some of 
these techniques were subsequently elaborated by Bogachev and Zarbaliev 
[6, 7] and also by Zarbaliev in his PhD thesis [40], however a complete proof 
has not been published as yet. 
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Remark 1.1. A polygonal line F £ IJn can be viewed as a vector sum of 
its consecutive edges, resulting in a given integer vector n = (ni, 712); due to 
the convexity property, the order of parts in the sum is uniquely determined. 
Hence, any such F represents an integer vector partition of n S Z^ or, more 
precisely, a strict vector partition (i.e., without proportional parts; see [32]). 
This observation incorporates the topic of convex lattice polygonal lines in a 
general theory of integer partitions. For ordinary, one-dimensional partitions, 
the problem of limit shape can also be set up, but for a special geometric 
object associated with partitions, called Young diagrams [33, 34]. 

1.2. Main result. Vershik [32, p. 20] pointed out that it would be inter- 
esting to study asymptotic properties of convex lattice polygonal lines under 
other probability measures Pn on i7„, and conjectured that the limit shape 
might be universal for some classes of measures. Independently, a similar 
hypothesis was put forward by Prokhorov [27]. 

In the present paper, we prove Vershik-Prokhorov's universality conjec- 
ture for a parametric family of probability measures P^ (0 < r < 00) on i7„ 
defined by 

(1-5) P;(r):=^, rG77„, 

with 

(1.6) b^\F):=llb[^, B::= J] b^\F), 

where the product is taken over all edges e^ of -T G 7T„, ki is the number of 
lattice points on the edge e^ except its left endpoint, and 

(1.7) ,j^^(r + .-l^^r(^+lV_(r+^-l)_ ^^^^^^^ 

Note that for r = 1 the measure (1.5) is reduced to the uniform distribu- 
tion on i7„. Qualitatively, formulas (1.6), (1.7) introduce certain probability 
weights for random edges on F by encouraging (r > 1) or discouraging 
(r < 1) lattice points on each edge as compared to the reference case r = 1. 
Assume that < ci < n2/ni < C2 < 00, and consider the standard 
scaling transformation Sn{x) = (x\ln\^X2ln2)^ x = (a;i,a;2) G K^. It is 
convenient to work with a sup-distance between the scaled polygonal lines 
Fn '■= Sn{F) {F G 7T„) and the limit curve 7*, based on the tangential 
parameterization of convex paths (see the Appendix, Section A.l). More 
specifically, for t G [0, 00] denote by in{t) the right endpoint of that part of 
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Fn where the tangent slope (wherever it exists) does not exceed t. Similarly, 
the tangential parameterization of the parabola arc 7* (see (1-3)) is given 
by the vector function 

The tangential distance between Fn and 7* is then defined as 
(1.9) drirn,l*)-= sup |e„(i)-<7*(t)|, 

0<t<oo 

where, as before, | • | is the Euclidean vector norm in M^ (cf. general definition 
(A. 3) of the metric d-i-{-, •) in the Appendix, Section A.l). 

We can now state our main result about the universality of the limit shape 
7* under the measures P^ (cf. Theorem 8.2). 

Theorem 1.1. For each r G (0,oo) and any e > 0, 

hm P„-{dr(/'n,7*)<4 = l- 
n— >oo 

It can be shown (see the Appendix, Section A.l) that the Hausdorff dis- 
tance d-^ (see (1.2)) is dominated by the tangential distance dj- defined in 
(A. 3) (however, these metrics are not equivalent). In particular. Theorem 1.1 
with r = 1 recovers the limit shape result (1.4) for the uniform distribution 
on Fin- As was mentioned above, in the original paper by Sinai [29] the proof 
of the limit shape result was only sketched, so even in the uniform case our 
proof seems to be the first complete implementation of Sinai's probabilistic 
method (which, as we will try to explain below, is far from straightforward). 

Let us also point out that Theorem 1.1 is a non-trivial extension of (1.4) 
since the measures P^ (r 7^ 1) are not close to the uniform distribution P^ 
in total variation distance (denoted by || • ||Ty)> and in fact ||P^ — P^HT^y ~^ 1 
as n — >• 00 (see Theorem A. 4 in the Appendix). 

The result of Theorem 1.1 for "pure" measures P^ readily extends to 
mixed measures. 

Theorem 1.2. Let p be a probability measure on (0,oo), and set 

/■oo 

(1.10) pan--= p^inpidr), Fen^. 

Jo 
Fhen, for any e > 0, 

hm pp{dr{rn,Y)<e} = i. 
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Proof. Follows from equation (1-10) and Theorem 1.1 by Lebesgue's 
dominated convergence theorem. D 

Theorem 1.2 shows that the limit shape result holds true (with the same 
limit 7*) when the parameter r specifying the distribution P^ is chosen 
at random. Using the terminology designed for settings with random en- 
vironments, Theorems 1.1 and 1.2 may be interpreted as "quenched" and 
"annealed" statements, respectively. 

Remark 1.2. The universality of the limit shape 7*, established in The- 
orem 1.1, is not a general rule but rather an exception, holding for some, 
but not all, probability measures on the polygonal space i7„. In fact, as was 
shown by Bogachev and Zarbaliev [7, 8], any C^-smooth, strictly convex 
curve 7 started at the origin may appear as the limit shape with respect to 
some probability measure Pn on 7T„, as n — ?• 00. 

Remark 1.3. The main results of the present paper have been recently 
reported (without proofs) in a brief note [9]. 

1.3. Methods. Our proof of Theorem 1.1 employs an elegant probabilis- 
tic approach first applied to convex lattice polygonal lines by Sinai [29]. 
This method is based on randomization of the (fixed) right endpoint ^ = n 
of polygonal lines P S 77^,, leading to the interpretation of a given (e.g., 
uniform) measure P„ on i7„ as the conditional distribution induced by a 
suitable probability measure Qz (depending on an auxiliary "free" parame- 
ter z = {zi,Z2)) defined on the "global" space 77 = Un77„ of all convex lattice 
polygonal lines (with finitely many edges). To make the measure Qz closer 
to Pn on the subspace i7„ C U specified by the condition ^ = n, it is natural 
to pick the parameter z from the asymptotic equation Ez{Cj = n{l + o(l)) 
(n —7- 00). Then, in principle, asymptotic properties of polygonal lines P 
(e.g., the limit shape) can be established first for {n,Qz) and then trans- 
ferred to (Un, Pn) via conditioning with respect to 77„ and using an appro- 
priate local limit theorem for the probability Qz{£, = n}. A great advantage 
of working with the measure Qz is that it may be chosen as a "multiplicative 
statistic" [34, 35] (i.e., a direct product of one-dimensional probability mea- 
sures), thus corresponding to the distribution of a sequence of independent 
random variables, which immediately brings in insights and well-developed 
analytical tools of probability theory. 

Sinai's approach in [29] was motivated by a heuristic analogy with sta- 
tistical mechanics, where similar ideas are well known in the context of 
asymptotic equivalence, in the thermodynamic limit, of various statistical 
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ensembles (i.e., microcanonical, canonical, and grand canonical) that may be 
associated with a given physical system (e.g., gas) by optional fixing of the 
total energy and/or the number of particles (see Ruelle [28]). In particular, 
Khinchin [22, 23] has pioneered a systematic use of local limit theorems of 
probability theory in problems of statistical mechanics. Deep connections 
between statistical properties of quantum systems (where discrete random 
structures naturally arise due to quantization) and asymptotic theory of ran- 
dom integer partitions are discussed in a series of papers by Vershik [34, 35] 
(see also the recent work by Comtet et al. [11] and further references therein). 
Note also that a general idea of randomization has proved instrumental in a 
large variety of combinatorial problems (see, e.g., [1, 2, 14, 15, 16, 24, 26, 33] 
and the vast bibliography therein). 

The probabilistic method is very insightful and efficient, as it makes the 
arguments heuristically transparent and natural. However, the practical im- 
plementation of this approach requires substantial work, especially in the 
two-dimensional context of convex lattice polygonal lines as compared to the 
one-dimensional case exemplified by integer partitions and the correspond- 
ing Young diagrams [33, 34]. To begin with, evaluation of expected values 
and some higher-order statistical moments of random polygonal lines leads 
one to deal with various sums over the set X of points x = (xi,a;2) G Z^ 
with co-prime coordinates (see Section 2.1). Sinai [29] was able to obtain 
the limit of some basic sums of such a kind by appealing to the known 
asymptotic density of the set X in Z^ (given by 6/7r^); however, this argu- 
mentation is insufficient for more refined asymptotics. In the present paper, 
we handle this technical problem by using the Mobius inversion formula (see 
Section 3), which enables one to reduce sums over X to more regular sums. 

As already mentioned, another crucial ingredient required for the prob- 
abilistic method is a suitable local limit theorem that furnishes a "bridge" 
between the global distribution Qz and the conditional one, P„. Analytical 
difficulties encountered in the proof of such a result are already significant 
in the case of ordinary integer partitions (for more details and concrete 
examples, see [2, 14, 13, 15, 16] and further references therein). The case 
of convex lattice polygonal lines, corresponding to two-dimensional strict 
vector partitions (see Remark 1.1), is notoriously tedious, even though the 
standard method of characteristic functions is still applicable. To the best 
of our knowledge, after the original paper by Sinai [29] where the result was 
just stated (with a minor error in the determinant of the covariance matrix 
[29, p. Ill]), full details have not been worked out in the literature (however, 
see [40]). We prove the following result in this direction (cf. Theorem 7.1). 
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Theorem 1.3. Suppose that the parameter z is chosen so that a^ 
SJ'(^) = n (1 + 0(1)). Then, as n^ 00, 



where K^ := Cov(^, ^) is the covariance matrix of the random vector^ {with 
respect to the probability measure Q^). 



Remark 1.4. The quantities a^ and K^, obtained via the measure Q^, 
depend in general on the parameter r as well. For the sake of notational 
convenience the latter is omitted, which should cause no confusion since r 
is always fixed, unless stated explicitly otherwise. 

One can show that the covariance matrix Kz is of the order of |np^, and 
in particular det K^ ~ const (72171-2)^'^ and \\Kz || = 0(|n|~^'^). From 
the right-hand side of (l-H), it is then clear that one needs to refine the 
error term in the asymptotic relation El{^) = n(l + o(l)) by estimating 
the deviation EK^) — n to at least the order of Inp'"^. We have been able to 
obtain the following estimate (cf. Theorem 5.1). 

Theorem 1.4. Under the conditions of Theorem 1.3, 

(1.12) E^(e) = n + o(|np/^), 71 ^00. 

The proof of this result is quite involved. The main idea is to apply the 
Mellin transform and use the inversion formula to obtain a suitable integral 
representation for the difference E'^^ — n oi the form (j = 1, 2) 

where Inzj happens to be of the order of |?i|~^'^ (according to the "optimal" 
choice of z as explained at the beginning of Section 1.3; cf. Theorem 3.1), 
Fj{s) is an explicit function analytic in the strip 1 < JRs < 2, and C(s) is the 
Riemann zeta function. As usual, to obtain a better estimate of the integral 
one has to shift the integration contour in (1.13) as far to the left as possible, 
and it turns out that to get an estimate of order o(|72p ^) one needs to enter 
the critical strip < Ks < 1, which requires information about zeroes of the 
zeta function in view of the denominator C,{s) in (1.13). 
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Layout. The rest of the paper is organized as follows. In Section 2, we 
explain the basics of the probability method in the polygonal context and 
define the parametric families of measures Q^ and P^ (0 < r < oo). In 
Section 3, we choose suitable values of the parameter z = (^1,22) (Theo- 
rem 3.1), which implies convergence of "expected" polygonal lines to the 
limit curve 7* (Section 4, Theorems 4.1 and 4.2). The refined error estimate 
(1.12) is proved in Section 5 (Theorem 5.1). Higher-order moment sums are 
analyzed in Section 6; in particular, the asymptotics of the covariance matrix 
Kz is obtained in Theorem 6.1. Section 7 is devoted to the proof of the local 
central limit theorem (Theorem 7.1). Finally, the limit shape result, with re- 
spect to both Ql and P^, is proved in Section 8 (Theorems 8.1 and 8.2). The 
Appendix includes necessary details about the tangential parameterization 
and the tangential metric dj- on the space of convex paths (Section A.l), 
as well as a discussion of the total variation distance between the measures 
P^ (r 7^ 1) and the uniform distribution P^ (Section A. 2, Theorems A. 2 
and A. 4). 

Notation. Let us fix some general notations frequently used in the paper. 
For a row-vector x = (3;i,2;2) € M^, its Euclidean norm (length) is denoted 
by |x| := {xf + Xg)^'^, and {x, y) := xy^ = xiyi + ^22/2 is the corresponding 
inner product of vectors x,y G M^. We denote Z+ := {A; G Z : fc > 0}, 

Z2_ := Z+ X Z+, and simialrly IR+ := {x G M : 2; > 0}, M^ := E+ x M+. 

2. Probability measures on spaces of convex polygonal lines. 

2.1. Encoding. As was observed by Sinai [29], one can encode convex lat- 
tice polygonal lines via suitable integer-valued functions. More specifically, 
consider the set X of all pairs of co-prime non- negative integers, 

(2.1) X:={x = (xi, X2) G 1?+ : gcd(xi, ^2) = 1}, 

where gcd(-,-) stands for the greatest common divisor of two integers. (In 
particular, the pairs (0,1) and (1,0) are included in this set, while (0,0) 
is not.) Let (p := (Z_|_)'^ be the space of functions on X with non-negative 
integer values, and consider the subspace of functions with finite support, 

<Po := {i^ e -P : #(supp z>) < 00}, 

where supp z/ := {x £ X : iy{x) > 0}. It is easy to see that the space ^0 
is in one-to-one correspondence with the space 77 = U^g^2 77.„ of all (finite) 
convex lattice polygonal lines: 

^o3 ly i — > Pen. 
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Indeed, let us interpret points x £ X as radius-vectors (pointing from the 
origin to x). Now, for any v £ (pQ, a finite collection of nonzero vectors 
{xv^x), X G supp i^}, arranged in the order of increase of their slope X2/X1 G 
[0,00], determines consecutive edges of some convex lattice polygonal line 
-T € 77. Conversely, vector edges of a lattice polygonal line F G 11 can be 
uniquely represented in the form xk, with x £ X and integer k > 0; setting 
i^(x) := k for such x and zero otherwise, we obtain a function v £ ^o- 
(The special case where i^ix) = for all x G Af corresponds to the "trivial" 
polygonal line Iq with coinciding endpoints.) 

That is to say, each x G X determines the direction of a potential edge, 
only utilized if x £ supp v, in which case the value z^(x) > specifies the 
scaling factor, altogether yielding a vector edge xv{x)] finally, assembling 
all such edges into a polygonal line is uniquely deterinined by the fixation 
of the starting point (at the origin) and the convexity property. 

Note that, according to the above construction, v{x) has the meaning of 
the number of lattice points on the edge xv{x) (except its left endpoint). 
The right endpoint ^ = ^r oi the polygonal line F £ 11 associated with a 
configuration u £ ^q\s expressed by the formula 

(2.2) e = 5^xKx). 

In what follows, we shall identify the spaces 77 and ^o- In particular, any 
probability measure on 77 can be treated as the distribution of a Z_|_-valued 
random field i^(-) on X with almost surely (a.s.) finite support. 

2.2. Global measure Qz and conditional measure Pn- Let bQ,bi,b2, ■ ■ ■ 
be a sequence of non-negative numbers such that 60 > (without loss of 
generality, we put bo = 1) and not all bk vanish for A; > 1, and assume that 
the generating function 

00 

(2.3) f3{s):=^bks' 

k=0 

is finite for |s| < 1. Let z = (21,^2) be a two-dimensional parameter, with 
zi,Z2 G (0, 1). Throughout the paper, we shall use the multi-index notation 

z^ := Zi^Z2^, X = (3;i,3;2) G Z^. 

We now define the "global" probability measure Qz on the space «? = 
(Z+) as the distribution of a random field z^ = {i^{x), x G X} with mutually 
independent values and marginal distributions of the form 

(2.4) Q^{j,(x) = k} = ^^, kez+ (xex). 

p(^^ ) 
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Proposition 2.1. For each z G (0, 1)^, the condition 
(2.5) m := n /5(^") < oo 

is necessary and sufficient in order that Qz{^o) = 1- 

Proof. According to (2.4), Qz{iy{x) > 0} = 1 - /3(z^)-^ (x G X). Since 
the random variables z^(x) are mutually independent for different x £ X, 
Borel-Cantelli's lemma implies that Qz{j-' G ^o} = 1 if and only if 

In turn, the latter inequality is equivalent to (2.5). D 

That is to say, under condition (2.5) a sample configuration of the ran- 
dom field i^(-) belongs (Q^-a.s.) to the space ^o and therefore determines a 
(random) finite polygonal line F £ IT. By the mutual independence of the 
values i^(x), the corresponding Q^-probability is given by 

where S, = Xl^eA' xv{x) is the right endpoint of F (see (2.2)), and 

(2.7) h{F) := l[ 6,(,) < oo, T G 77. 

xex 

Note that the infinite product in (2.7) contains only finitely many terms dif- 
ferent from 1, since for x ^ supp ly we have b,yi^x) = ^o = 1- Hence, expression 

(2.7) can be rewritten in a more intrinsic form (cf. (1.6)) 

(2.8) b{F) := n h,, 

where the product is taken over all edges Cj of 7^ G 77„ and ki is the number 
of lattice points on the edge e, except its left endpoint (see Section 2.1). 

In particular, for the trivial polygonal line 7o -H- z^ = formula (2.6) 
yields 

Q,(ro) = ^(z)-^>o. 

Note, however, that Qz{Fq) < 1 since, due to our assumptions, (2.3) implies 
/3(s) > /3(0) = 1 for s > and hence, according to (2.5), /3{z) > 1. 
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On the subspace 77„, C 11 of polygonal lines with the right endpoint fixed 
at n = (711,722), the measure Qz induces the conditional distribution 

(2.9) p^{r):=Qz{r\nn) = §^^, r e n^, 

provided, of course, that QziHn) > (i.e., there is at least one F G i7.„ with 
b{r) > 0, cf. (2.6)). The parameter z may be dropped from the notation for 
Pn due to the following fact. 

Proposition 2.2. The measure Pn in (2.9) does not depend on z. 

Proof. U Hn B P -^ i^ £ ^o then ^ = n and hence formula (2.6) is 
reduced to 

Q.{P) = m^, PGPn. 
/3{z) 

Accordingly, using (2.5) and (2.9) we get the expression 
(2.10) Pn{P) = ^ ^^^[,^,. , r G 7T„, 

which is z-free. D 

2.3. Parametric families {Q^} and {P^}- Let us consider a special para- 
metric family of measures {Q^, < r < oo}, determined by formula (2.4) 
with the coefficients 6^ of the form 

(2.11) ,.^=(^ + ^-i) = !l!:±i)_i!:±i^, ,,^^ 

(note that 6g = (^^ q ) = 1, in accordance with our convention in Section 2.2). 
By the binomial expansion formula, the generating function (2.3) of the 
sequence (2.11) is given by 

(2.12) /5'■(s) = (l-.)-^ |.|<1, 

and from (2.4) it follows that under the law Q^ the random variable z^(x) 
has the probability generating function 

(2.13) j;r(,K^h^^(f£!l^ (1-^T o<s<l. 
Consequently, formula (2.4) specializes to 

(2.14) Q:4z.(x) = A;} = r^'^^" ^"j 2^=^(1 -z^)^ k e Z+ {x € X). 
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That is to say, with respect to the measure Q^ the random variable i^(x) 
has a negative binomial distribution with parameters r and p = 1 — z^ [12, 
§VI.8, p. 165]; in particular, its expected value and variance are given by 
(see [12, §XI.2, p. 269]) 

(2.15) ElHx)] = -^, Var[z.(x)] = ^^^^ (x G X). 

According to formulas (2.9) and (2.10), the corresponding conditional 
measure -P^(-) := (5^(-|-/7„) is expressed as 

where V{r) is given by the general formula (2.8) specialized to the coeffi- 
cients 6^ defined in (2.11). 

In the special case r = 1, we have h\ = (fc) = 1 so that (2.14) is reduced 
to the geometric distribution (with parameter p = 1 — z^) 

QlMx) = k} = z^'^il - z^), keZ+ (x G X), 

whereas the conditional measure (2.16) specifies the uniform distribution on 
Hn (cf. [29]): 

^»(^> = #(ky ''-"- 

Remark 2.1. Since 6^+1/6^ = (r + k)/{k + 1), the sequence {6^} is 
strictly increasing or decreasing in k according as r > 1 or r < 1, respectively. 
That is to say, the measures Q^ ^^d P^ encourage (r > 1) or discourage 
(r < 1) lattice points on edges, as compared to the reference case r = 1. 

It is easy to see that condition (2.5) is satisfied and, by Proposition 2.1, 
Qli'Po) = 1, < r < 00. 

Indeed, using (2.12) we have 



p'{z) = n (1 - ^")"' = e^pf -^ E 1^(1 - ^") ) 



< 00 



whenever "^^f^x lii(l ~ ^^) > —00, and the latter condition is fulfilled since 

00 00 ^ 

x&X ^&2 ^,=0 ^^=0 Vi ^IJU Z2) 
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3. Calibration of the parameter z. In what follows, the asymptotic 
notation of the form a„ x 6„ (where n = (ni, 112)) means that 

< liminf — ^ < limsup -— < 00. 

We also use the standard notation a„ ~ 6„, for an/bn — )• 1 as ni,n2 — )■ 00. 

Throughout the paper, we shall work under the following convention 
about the limit n = (ni,?i2) ~^ 00. 

Assumption 3.1. The notation n — )■ oo signifies that ni,n2 — )• oo in 
such a way that ni x n2. In particular, this implies that \n\ = {nl + n^Y''^ — t- 
oo as n — 7- oo, and ni x \n\, 722 x |n|. 

The goal of this section is to use the freedom of the conditional distribu- 
tion Pn{-) = Q^{'\nn) from the parameter z (see Proposition 2.2) in order 
to better adapt the measure Q^ to the subspace 77„ C U determined by the 
condition ^ = n (where ^ = (^11^2) is defined in (2.2)). To this end, it is 
natural to require that the latter condition be satisfied (at least asymptot- 
ically) for the expected value of S, (cf. [29, 6]). More precisely, we will seek 
z = {zi,Z2) as a solution to the following asymptotic equations: 

(3.1) i^r(6)~ni, El{i2)^n2 (n^oo), 

where E'J denotes expectation with respect to the distribution Q^- 
From (2.2), using the first formula in (2.15), we obtain 

(3.2) ^r(e) = E-T^ = ^EE--'^- 

x£X k=lxeX 

Let us represent the parameters zi, Z2 in the form 

(3.3) ^i = e-°^ aj = 6jnj^/^ (i = l,2), 

where the quantities 61,62 > (possibly depending on the ratio n2/ni) are 
presumed to be bounded from above and separated from zero. Hence, (3.2) 
takes the form 

00 

k{a,x) 



(3.4) ^1(6 = ^EE 



xe 

k=lx€X 
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Theorem 3.1. Conditions (3.1) are satisfied if Si, 62 in (3.3) are chosen 
to be 

(3.5) 61 = Kr^'\n2/nif/\ 82 = Kr^''\ni/n2f'\ 

where k := (^(3)/C(2)) and C{s) = Y1T=1^~'^ ^•^ ^^^ Riemann zeta func- 
tion. 

Proof. Let us prove the first of the asymptotic relations (3.1). Set 

(3.6) /(x) :=rxie-<"'^\ a; G M^, 

and 

F\h) := Y, fihx), h>0. 

Then we can rewrite (3.4) in projection to the first coordinate as 

k=l x&X k=l 

Let us also consider the function 

00 00 

(3.8) F{h) := Y^ F^hm) =YY1 /(^"^^)' ^ > 

m=l m=0 x£X 

(adding terms with m = does not affect the sum, since /(•) vanishes at 
the origin). Recalling the definition of the set X (see (2.1)), we note that Z^ 
can be decomposed as a disjoint union of multiples of X: Z^ = [JJ^^q "^'^• 
Hence, the double sum in (3.8) is reduced to 

00 00 

(3.9) xeZ^ a:i = l X2=0 



rhe 



-ha 



By the Mobius inversion formula (see [18, Theorem 270, p. 237]), 

00 00 

(3.10) F{h) = Y F^{hm) ^^ F\h) = ^ /i(m)F(/im), 

m=l m=l 

where /i(?n,) {m G N) is the Mobius function defined as follows: /u(l) = 1, 
IJ,{m) = (— 1)*^ if m is a product of d different primes, and iJ.{m) = if m 
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has a squared factor [18, §16.3, p. 234]; in particular, |/u(-)| < 1. A sufficient 
condition for (3.10) is that the double series Xlfc m l-^''(^^"'')l should be 
convergent, which is easily verified in our case: F*[-) > and, according to 
(3.8) and (3.9), 

oo CO oo ^ —hk(y-\ 

k,m=l k=l k=l ^ ^ ^ ' 

Using (3.9) and (3.10), we can rewrite (3.7) as 
(3.11) 



OO ^ OO 



k=\ m=l k,m=l ^ ^ ^ ' 

Note that (3.3) and (3.5) imply 

2 ^'^^ 2 ''^^ 

(3.12) aiQ2 = , aiUo = , a2n2 = aini, 

ni 71-2 

where k is defined in Theorem 3.1. Hence, we can rewrite (3.11) in the form 

(3.13) n-^EZii,) = \y rnMm)afa2e-^-^ 

k,m=l ^ ^ ^ ' 

We now need an elementary estimate, which will also be instrumental 
later on. 

Lemma 3.2. For any k > 0, 9 > 0, there exists C = C{k,0) > such 
that, for all t > 0, 

(3-14) T^ TTT < ^^ • 

y I (1 - e-*)'= - t^ 

Proof of Lemma 3.2. Set g{t) := t'^e~^*/^(l - e~*)~'^ and note that 
lim g{t) = 1, lim g{t) = 0. 

t— !>0+ t— >oo 

By continuity, the function g{t) is bounded on (0, oo), and (3.14) follows. D 

By Lemma 3.2, the general term of the series (3.13) is estimated, uniformly 
in k and m, by 0[k~^m~^). Hence, by Lebesgue's dominated convergence 
theorem one can pass to the limit in (3.13) termwise: 



^ oo ^ oo 

(3.15) li^n^^E:{Ci) = -J2vsT. 

71— >oo K ■^ — ' K ^ — ' 



k=l m=\ ^^ ' 
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Here the expression for the second sum (over m) is obtained using the 
Mobius inversion formula (3.10) with F^{h) = /i"^ F{h) = Em=i(^"^)"^ = 
/i-2C(2) (cf. [18, Theorem 287, p. 250]). 
Similarly, we can check that, as n — )• oo, 

~ipr(f:^^J_ Y^ m/x(m)aia|e-^"""2 

k,m=l ^ ^^ ' 

The theorem is proved. D 

Remark 3.1. The term 1/C(2) = O/vr^ appearing in formula (3.15) and 
the like, equals the asymptotic density of co-prime pairs x = (xi,X2) G X 
among all integer points on Z^ (see [18, Theorem 459, p. 409]). 

Assumption 3.2. Throughout the rest of the paper, we assume that the 
parameters zi, Z2 are chosen according to formulas (3.3), (3.5). In particular, 
the measure Q^ becomes dependent on n = (77,1,77-2), as well as all Q^- 
probabilities and mean values. 

4. Asymptotics of "expected" polygonal lines. For F £ 11, denote 
by F{t) {t G [0,00]) the part of F where the slope does not exceed tn2/ni. 
Hence, the path Fn{t) = Sn{F{t)) serves as a tangential parameterization 
of the scaled polygonal line Fn = Sn{F), where Snixi,X2) = {xi/ni,X2/n2) 
(see Section 1.2, and also Section A.l below). Consider the set 

(4.1) X{t) ■.= {x£ X : X2IXX < tn2/ni}, t G [0, 00]. 

According to the association TI 3 F ■v^ u £ <^q described in Section 2.1, 
for each t G [0, cxo] the polygonal line F{t) is determined by a truncated 
configuration {^{x), x G X{t)}, hence its right endpoint ^(t) = {ii{t) , £,2{f)) 
is given by 

(4.2) m= Yl ^^(^)' ie[0,oo]. 

x£X(t) 

In particular, A'(oo) = X, ^(00) = i (cf. (2.2)). Similarly to (3.2) and (3.4), 



00 



(4.3) ElUt)] =rY.H ^e-^<"'-), t G [0,oo]. 

k=lx£X{t) 

Let us also set (cf. (1.8)) 

(4.4) 9Kt):=:^±^, ^2W:=(Y^, tG[0,oo]. 
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As will be verified in the Appendix (see Section A.l), the vector- function 
g*{t) = (fi'i(i),52(^)) gives a tangential parameterization of the parabola 7* 
defined in (1.3). 

The goal of this section is to establish the convergence of the (scaled) 
expectation £'^[^(t)] to the limit g*{t), first for each t G [0,oo] (Section 4.1) 
and then uniformly in i G [0,oo] (Section 4.2). 

4.1. Pointwise convergence. 

Theorem 4.1. For each t G [0, oo], 
(4.5) hm nj'Eli^jit)] = g*{t) {j = 1, 2). 



Proof. Theorem 3.1 implies that (4.5) holds for t = oo. Assume that 
t < oo and let j = 1 (the case j = 2 is considered in a similar manner). 
Setting for brevity c„ := n^jnx and arguing as in the proof of Theorem 3.1 
(see (3.4), (3.7) and (3.13)), from (4.3) we obtain 



^2 



k,m=l xi=l 2:2=0 

00 / ^ 00 

(4.6) = r Y "^ Trig, E ^1^'""""'"' (1 - e-'='""^(^^+i)' 

k,m=l xi=l 

where £2 = X2{t) denotes the integer part of tc^xi, so that 

(4.7) < tc„xi -X2 < 1. 

Aiming to replace X2 + 1 by tc,«xi in (4.6), we recall (3.12) and rewrite the 
sum over xi as 

00 

(4.8) Y ^1^'""""'"' (1 - e-*^™"i*^i) + Afc,^(t, a), 

where 

00 

Afc,™(t,a) := Y Xie-''""''''^^^+^Ul - Q-kma2{x2+l-tc„Xi)\ ^ 

Xl=l 

Using that < X2 + 1 — ic„xi < 1 (see (4.7)) and applying Lemma 3.2, we 
obtain, uniformly in A;,rn- > 1 and t G [0, 00], 

A, (f n) ^ , p-kmai -mai/2 

< f^''"7"^ < y xie-'^"^-!^^ = -^—^ -. = 0(1) ^ ^ . 

Xl = l 



18 L. V. BOGACHEV AND S. M. ZARBALIEV 

Substituting this estimate into (4.6), we see that the error resulting from 
the replacement of £2 + 1 by ic„xi is dominated by 

OK^)E^E^^;^ = OK^)ln(l-e--/^)=0(ar^lnaO. 

fc=l m=l 

Returning to representation (4.6) and computing the sum in (4.8), we find 



(4.9) E:Mt)] = ry '"^T • " /\, 

V / 2Lsx\ n / J -^ _ a—kma.2 (\ — p—kmoLiy\2 



k,m=l 



y=l+t 



Passing to the limit by Lebesgue's dominated convergence theorem, similarly 
to the proof of Theorem 3.1 (cf. (3.15)) we get, as n — )• 00, 



( 1 m — 1 \ \ ' / / 



t^ + 2t 



2 ' 



which coincides with gKt), as claimed. D 

4.2. Uniform convergence. There is a stronger version of Theorem 4.1. 

Theorem 4.2. Convergence in (4.5) is uniform in t £ [0,oo], that is, 
lim sup I nfEl [Cj (t)] - g* (t) | = (i = 1, 2) . 

For the proof, we need the following general lemma. 

Lemma 4.3. Let {fnit)} be a sequence of nondecreasing functions on 
a finite interval [a,b], such that, for each t € [a,b], hnin-^oo fn{t) = fit), 
where f{t) is a continuous (nondecreasing) function on [a,b]. Then the 
convergence fn{t) — )• f{t) as n ^ 00 is uniform on [a,b]. 

Proof of Lemma 4.3. Since / is continuous on a closed interval [a, b], it 
is uniformly continuous. Therefore, for any e > there exists 5 > such that 
\f{t') - f{t)\ < e whenever \t' -t\< 5. Let a = to < *i < • • • < % = ^ be a 
partition such that maxi<j<7v (ti — tj-i) < 5. Since lim„_!.oo fn{ti) = f{ti) for 
each i = Q,1, . . . ,N , there exists n* such that maxo<i<Ar \fn{ti) — f{ti)\ < £ 
for all n > n* . By monotonicity of fn and /, this implies that for any 
t G [ij,tj+i] and all n > n* 

fnit) - fit) < fniti+i) - fiti) < /„(t.+i) - fiti+i) + e<2e. 
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Similarly, /„,(t) - f{t) > -2e. Therefore, sup^gj^;,] |/„(t) - f{t)\ < 2e, and 
the uniform convergence follows. D 

Proof of Theorem 4.2. Suppose that j = I (the case j = 2 is handled 
similarly). Note that for each n the function 

Til ^-^ 

x€X(t) 

is nondecreasing in t. Therefore, by Lemma 4.3 the convergence (4.5) is 
uniform on any interval [0,t*]. Furthermore, since n^ El[^i{oo)] — t- gKoo) 
and the function gl{t) is continuous at infinity (see (4.4)), for the proof of 
the uniform convergence on a suitable interval [t* , oo] it suffices to show that 
for any e > one can choose t* such that, for all large enough ni, n2 and 
aU t > t*, 

(4.10) n^'El\Ci{oo)-^i{t)\<e. 
On account of (4.9) we have 

(4.11) ElUoo)-m]= ±-P^ ^"'"^"^^''"' , 

k,m=l K^ ^ ) 

+ 0{a^^ In ai). 
Note that by Lemma 3.2, uniformly m k,m> 1, 

g-fcm«2 g-fcmai(l+i) 0(1) 

1 - e-'^'^^a ' (i_e-fc"^°i(i+*))2 ~ a^a2(A:m)3(l + t)2' 
Returning to (4.11), we obtain, uniformly in t > t*, 

2 ^rr. . N . . M 0(1) ^^ 1 ^22, 1 0(1) 

'^ ' ^ k=l m=l ^ ' 

whence by (3.3) we get (4.10). D 

5. Further refinement. For future applications, we need to refine the 
asymptotic formulas (3.1) by estimating the error term. The following the- 
orem is one of the main technical ingredients of our work. 

Theorem 5.1. Suppose that the parameter z is chosen according to for- 
mulas (3.3), (3.5), so that El{^) = n(l + o(l)) (see Theorem 3.1). Then 
^ziC) = n + o(|?ip/^) as n -^ oo. 

For the proof of this theorem, some preparations are needed. 
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5.1. Approximation of sums by integrals. Let a function /: M^ — > M 
be continuous and absolutely integrable on M^, together with its partial 
derivatives up to the second order. Set 

(5.1) Fih) := Y^ f{hx), h > 

(as verified below, the series in (5.1) is absolutely convergent for all h > 0), 
and assume that for some /3 > 2 

(5.2) F{h) = 0{h-P), h^oo. 

Consider the Mellin transform of F{h) (see, e.g., [38, Ch. VI, §9]), 

POO 

(5.3) F{s) := / h'-^F{h)dh. 

Jo 

Lemma 5.2. Under the above conditions, the function F(s) is meromor- 
phic in the strip 1 < Ks < /?, with a single {simple) pole at s = 2. Moreover, 
F{s) satisfies the identity 

POO 

(5.4) Fis)= h'^^Af{h)dh, l<^s<2, 

Jo 

where 

(5.5) Af{h) := F{h) - ^ [ f{x)dx, h > 0. 

"- Jr^_^ 

Remark 5.1. Identity (5.4) is a two-dimensional analogue of Miintz's 
formula for univariate functions (see [31, §2.11, pp. 28-29]). 

Proof of Lemma 5.2. Let a function (p : M+ — ^ R be continuous and 
continuously differentiable, and suppose that both (p and cf)' are absolutely 
integrable on M+. It follows that lim2,.^oo 0(2;) = 0; indeed, note that 



POO PX 

/ (j)'{x)dx= lim / (j)'{y)dy= lim 
Jo ^^°° Jo ^'^"^ 



(x) - 0(0), 



hence lim^^^oo <Pix) exists and, since (j) is integrable, the limit must equal 
zero. Then the well-known Euler-Maclaurin summation formula states that 

00 -J /■ 00 /"OO 

(5.6) Yl <^(^-?') = h "^^^^ "^"^ + / ^^ (D "^'^"^^ '^'^' 

■ ri "/ U "^ U 
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where i?i(x) := x - [x] - 1 (cf. [5, §A.4, p. 254]). 

Applying formula (5.6) twice to the double series (5.1), we obtain 

f ~ / ^1 \ ~ / xo\ d f(x^ 
J^2 \ h J \ h J 3xi9x2 



Since |-Bi(-)| ^ !> the above conditions on the function / imply that all 
integrals in (5.7) exist, hence F{h) is well defined for all h > 0. Moreover, 
from (5.7) it follows that 

(5.8) F{h)=0{h~^), Af{h) = 0{h^^) (h^O). 

The estimates (5.2) and (5.8) imply that F{s) as defined in (5.3) is a 
regular function for 2 < 3f?s < /3. Let us now note that for such s we can 
rewrite (5.3) as 

h'^^F{h)dh+ / h'~^F{h)dh 

= / h'~^F{h)dh+ / h'-^dh / f{x)dx+ / h'~^Af{h)dh 
Ji Jo Jw^ Jo 

(5.9) = [ h'-^F(h)dh+^— I f{x)dx+ I h'^^Af(h)dh. 

J I s-2 J^2^ Jo 

According to condition (5.2), the first term on the right-hand side of (5.9), 
as a function of s, is regular for Ks < /3, whereas the last term is regular for 
Ks > 1 by (5.8). Hence, formula (5.9) furnishes an analytic continuation of 
the function F{s) into the strip 1 < Ks < /3, where it is meromorphic and, 
moreover, has a single (simple) pole at point s = 2. Finally, observing that 



oo 



^ ' h'~'^dh, Ks<2, 



s-2 ji 

and rearranging the terms in (5.9) using (5.5), we obtain (5.4). D 

Lemma 5.3. Under the conditions of Lemma 5.2, 
(5.10) Af{h) = -— h~'F{s)ds, 1< c < 2. 
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Proof. From (5.5), (5.7) we have Af{h) = 0{h^'^) as /i — ;> oo. Combined 
with the estimate (5.8) estabhshed in the proof of Lemma 5.2, this imphes 
that the integral (5.4) converges absolutely in the strip 1 < Ks < 2. The 
representation (5.10) then follows from (5.4) by the inversion formula for 
the Melhn transform (see [38, Theorem 9a, pp. 246-247]). D 

5.2. Proof of Theorem 5.1. Let us consider the first coordinate, ^i (for 
^2 the proof is similar). The proof consists of several steps. 

Step 1. According to (3.11) we have 



(5.11) E:{i,)= Y: ^Fikm), 

k,m=l 



fj,{m 
k,m=l 

where (see (3.6), (3.9)) 

(h'r\ = . 

(1 -e-°i^)2(l-e-"2'^) 



IT' hi (zt fJ^ 1 ' t 

F{h) = Y. f{h.) = rr—z^^^K^r—z^y h > 0, 



xgZ; 



2 



/(x) = r3;ie-<"'^\ xGR+. 



Note that 



j-OO l-OO 

f{x)dx = r xie-"^^'dxi e'^^^Mxa = ^— . 
Jo Jo ai"2 

Moreover, using (3.12) we have 

(cf. (3.15)). Subtracting (5.12) from (5.11), we obtain the representation 

oo 

(5.13) El{^i)-ni= E rnn{m) Af{km), 

k,m=l 

where Aj(/i) is defined in (5.5). Clearly, the functions / and F satisfy the 
hypotheses of Lemma 5.2 (with /3 = oo). Setting c„ := n2/ni and using 
(3.12), the Mellin transform of F{h) defined by (5.3) can be represented as 

Fis)=ra-^'+^^F{s), 
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where 

(^■") H^)-y„ (i_e4(l-e-./^-) '"'- '''*'■ 

As a result, applying Lemma 5.3 we can rewrite (5.13) as 
(5.15) 



fc,m=l 



S'tep 2. It is not difficult to find explicitly the analytic continuation of 
the function F[s) into the domain 1 < 3f?s < 2. To this end, let us represent 
the integral (5.14) as 

(5.16) Fis) = Jis)+c^j^ _l__d,+ _y^ (T^l!'. 
where 

(5.17) .Us) := r J^iLi. { ' ^ _£._!) d,. 

Jo (l-e-!')2 Vl-e-»/'. !/ 2/ 

The last two integrals in (5.16) are easily evaluated: 

r-oo ,,s— Ip— J/ /"OO °° °° ^oo 



Jo (i-e S'j Jo ;-^i ;-^i ./O 

°° 1 /"CXI 

(5.18) =Y.Ti^ u^-^e--du = as-i)ns). 

k=l ^ -^0 
where r(s) = f^ n'*~^e~" du is the gamma function, and similarly 

(5.19) y^ ^^^^dy = c(.)r(« + i). 

Substituting expressions (5.18) and (5.19) into (5.16), we obtain 
(5.20) F{s) = J{s) + CnCis - l)r(s) + ^ ({s)T{s + 1). 

Since the expression in the parentheses in (5.17) is 0{y) as y — )■ and 
0(1) as 2/ — )■ oo, the integral (5.17) is absolutely convergent (and therefore 
the function J{s) is regular) for ?R.s > 0. Furthermore, it is well known that 
r(s) is analytic for Ks > [30, §4.41, p. 148], while C(s) has a single pole 
at point s = 1 [30, §4.43, p. 152]. Hence, the right-hand side of (5.20) is 
meromorphic in the semi-plane 3?s > with poles at s = 1 and s = 2. 
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Step 3. Let us estimate the function F{c + it) as t — t- oo. Firstly, by 
integration by parts in (5.17) it is easy to show that, uniformly in a strip 
< ci < a < C2 < oo, 

(5.21) J{a + it) = 0{\t\~'^), t^oo. 

The gamma function in such a strip is known to satisfy a uniform estimate 

(5.22) r(a + ii) = 0(l)|ir-(^/2)g--|t|/2^ ^^^ 

(see [30, §4.42, p. 151]). Furthermore, the zeta function is obviously bounded 
in any semi-plane cr > ci > 1 : 

oo ^ oo ^ oo ^ 

(5.23) \Ci- + ^t)\<Eu^^\=^^-^^^= ^^'^- 

n=l n=l n=l 

We also have the following bounds, uniform in a, on the growth of the zeta 
function as t — )■ oo (see [20, Theorem 1.9, p. 25]): 



(5.24) ({a + it) = < 



f 0(ln|i|), 1 < C7 < 2, 

0[t^^-''y^ln\t\), 0<a<l, 
^0{t^^^~''ln\t\), a<0. 



As a result, by (5.22), (5.23) and (5.24) the second and third summands on 
the right-hand side of (5.20) give only exponentially small contributions as 
compared to (5.21), so that 

(5.25) F{c + it) = 0{\t\-'^), t^oo (1< c < 2). 

Step 4. In view of (5.25), for 1 < c < 2 there is an absolute convergence 
on the right-hand side of (5.15), 

y m|/i(m)| / ' ^^' ^ |dg| 



fc,m=l 



oo 



oo 



^^cTtT.T^iY.—cJ \F{c + ^t)\<it <oo. 

"l fc=l m=l "^ •^-°° 

Hence, the summation and integration in (5.15) can be interchanged to yield 



-'«"— i/rsss^E 



■ds 



IC-lOO "1 f.^^ ■" „^i "1 



(5.26) ^^r^-F(.K(»+l) 
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While evaluating the sum over m here, we used the Mobius inversion formula 
(3.10) with F»(/i) = /^-^ F{h) = Y^mi^^y = h~%{s) (cf. (3.15); see also 
[18, Theorem 287, p. 250]). Substituting (5.20) into (5.26), we finally obtain 



(5.27) 
where 

(5.28) 



Kiii) -^^i = lj- "^(s) ds (1< c < 2), 



^{s) :-- 



C(s + i) 



a 



s+l 



J{s)+Cnas-1)T{S) , 1^, ^^. 

0^^ + 2^(^ + 1) 



and the function J{s) is given by (5.17). 

Step 5. By the La Vallee Poussin theorem (see [21, §4.2, Theorem 5, 
p. 69]), there exists a constant A > such that ("(o" + it) 7^ in the domain 



(5.29) 



a> 1 



A 



In t +2 



=: r?(t), 



t e 



Moreover, it is known (see [31, Eq. (3.11.8), p. 60]) that in the domain (5.29) 
the following uniform estimate holds: 



(5.30) 



1 



0(ln|t|), 



t — )• 00. 



C{a + it) 
Without loss of generality, one can assume ^4 < ln2, so that (see (5.29)) 



r/(t) > ri{Q) = 1 - _- > 0, t G M. 

mi 

Therefore, ^{s) (see (5.28)) is regular for all s = a+it such that 2 > a > rj{t) 
{t G M). 

Let us show that the integration contour JRs = c in (5.27) can be replaced 
by the curve a = rj[t) {t G M). By the Cauchy theorem, it suffices to check 
that 

i-c+iT 

lim / l2/(s) ds = 0. 



We have 



(5.31) 



c+iT 



l^(s)ds 



r){T)+iT 



pc pc 

< \1'{a + iT)\da< {I'ia + iT)\da. 

JrjiT) JrjiO) 



26 



L. V. BOGACHEV AND S. M. ZARBALIEV 



In view of the remark after formula (5.30), we have r]{0) > 0, hence apph- 
cation of the estimate (5.23) gives, for s = a + iT, rj{T) < cr < c, 



C(5 + i) 



a 



s+l 



< C(^ + i) < C(^(o) + 1) 



a 



(T+l 



a 



c+l 



(since ai — )■ 0, we may assume that ai < 1). 

To estimate the expression in the square brackets in (5.28), we use the 
estimates (5.21), (5.22), (5.24) and (5.30). As a result, we obtain 



(5.32) 



^{a + iT) = 0{\T\~'^ln\T\), T ^ ±oo, 



which implies that the right-hand side of (5.31) tends to zero as T — )• ±00, 
as required. Therefore, the integral in (5.27) can be rewritten in the form 



(5.33) 



Dn-.-- 



<lr{r]{t) + it)d{i]{t) + it). 



Step 6. It remains to estimate the integral (5.33) as n — > cc. Let us set 
Ms) := af+^a/(s) 

= C(^ + i) 



(5.34) 



as) +2^(^+1) 



(see (5.28)), then equation (5.33) is rewritten as 

/oo 
aJ-''(*)-**lZ^o(r?(t) + it)(r?'(t) + i)di. 
-00 

Using that ai = 5i/n^'^ (see (3.3)), we get 

/oo 
al~'^^'^\%{v{t)+it)\{Wit)\ + l)dt 
-00 

/oo 
«!-''(*) |<^o(r?(i)+ii)|dt, 
-00 

since by (5.29) 



\v'{t)\ 



A 



< 



A 



{\t\ +2)ln^(|t| +2) ~ 21n^2 



0(1). 



Let us now note that, as n — ?• 00, the integrand function in (5.35) tends to 
zero for each t, because ai — )■ and 1 — ri{t) > (see (5.29)). Finally, eligi- 
bility of passing to the limit under the integral sign follows from Lebesgue's 
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dominated convergence theorem. Indeed, the integrand function in (5.35) is 
bounded by \^Q{r]{t) + it)\, and integrabihty of the latter is easily checked 
by applying the estimates (5.21), (5.22), (5.24) and (5.30) to the expression 
(5.34), which yields (cf. (5.32)) 

\%{7]{t) + it)\ =0{\t\-^ln\t\), t^±oo. 

Thus, we have shown that the integral in (5.35) is o(l) as n — )• oo, hence 
Dn = o(|np'^). Substituting this estimate into (5.27), we obtain the state- 
ment of Theorem 5.1. The proof is complete. 

6. Asymptotics of higher-order moments. 

6.1. The variance. According to the second formula in (2.15), we have 

(6.1) Var[Kx)] = .. ^ x)2 = ^E ^^'"- 

^ ^ ) k=i 

Let Kz := Cov(i^,i^) be the covariance matrix (with respect to the measure 
(5^) of the random vector ^ = '^^ex ■^'^i-'^)- Recalling that the random 
variables I'ix) are independent for different x £ X and using (6.1), we see 
that the elements Kz{i,j) = Cov{S,i,S,j) of the matrix K^ are given by 

oo 

(6.2) K,{i,j) = Y^ XiXjVar[iy{x)] = rJ2Yl ^^i^j^'"'^ ^^J ^ {1'2}- 

x&X k=lx£X 

Theorem 6.1. Asn^co, 

(6.3) Kz{i,j)^ ^''''^f^\ ,„ i,je{l,2}, 

where k is defined in Theorem 3.1 and the matrix B := (-By) is given by 
(0.4) sJ^n,/n, 1 Y 

Proof. Let us consider Kz{l, 1) (the other elements of K^ are analyzed 
in a similar manner). Substituting (3.3) into (6.2), we obtain 

oo 

(6.5) K,(l,l)=r^^fex?e-'=<°'^>. 

k=lx&X 
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Using the Mobius inversion formula (3.10), similarly to (3.13) expression 
(6.5) can be rewritten in the form 

oo 

i^^(l,l)=r Y^ km'^fi{m) Y, x?e-*^'"<"'^> 

k,rn=l xGZI 

oo oo oo 

kma2X2 



= r Y, km^Km) Y 2;?e-^™"i^i Y ^ 

k,m=l 2:1=1 X2=0 

(6.6) ^^y^nPj^- 



^2^-kma,x,_ 



k,m=l a;i=l 

Note also that 

;i + e-^™"i) 0(1) 



(6-7) Y 



2 —kmaixi 



Xie 



1 M — p—kmai\3 



a;i=l ^ 

Returning to the representation (6.6) and using (6.7), we obtain 

1,3 „, ^—kmaifi i^ „—kma 



kma2\ 

fc, 771=1 



(6.8) aia,K,{l,l)=r Y ^^"^'^-') ^i_\-k^..f^i_ ,- 



By Lemma 3.2, the general term in the series (6.8) admits a uniform estimate 
0(A;~^?n~^). Hence, by Lebesgue's dominated convergence theorem one can 
pass to the limit in (6.8) to obtain 

af 02^^(1,1) -^ ,' =2rK^ ai,a2^0. 

Using (3.3) and (3.5), this yields 

i^,(l,l)~^^^(nln2)'/^ n^oo, 

as required (cf. (6.3), (6.4)). D 

6.2. Statistical moments of u[x). Denote 

(6.9) vq{x) := v{x) - El[v{x)], x^X, 
and for A; G N set 

(6.10) mk{x) := El[v{x)% /ifc(x) := El\uo{xf\ 

(for notational simplicity, we suppress the dependence on r and z). 
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Lemma 6.2. For each k £N and all x £ X, 

(6.11) /ifc(x) < 2'=mfc(x). 

Proof. Omitting for brevity the argument x, by Newton's binomial for- 
mula and Lyapunov's inequality we obtain 






j=0 

and (6.11) follows. D 

Lemma 6.3. For each k £ N, there exist positive constants c^ = Ck{r) 
and Ck = Ck{r) such that, for all x £ X, 

(6.12) -^h^<m,ix)< ^"'^ 



(1-Z^)^- - «^ ^ - (l_2^)fc- 

Proof. Fix x £ JV and let ip{s) = (Pi,(x){s) ■= E"^ [e*'''^^^)] be the char- 
acteristic function of the random variable z^(x) with respect to the measure 
Ql. From (2.13) it follows that 

(0.13) ^is)-^^'""''- C-^")' 



/?^(z^) (1 - z^e'-'^y ■ 

Let us first prove that for any k £N 

(6 14) n-zn-^^-z^Vc-.^^^^^^L 

where Cj^^ = Cj^kix) > 0. Indeed, if A; = 1 then differentiation of (6.13) yields 



(l-.-)-^ 



ds (1 - z^e*'')''+i ' 

which is in accordance with (6.14) if we put ci,i := r. Assume now that 
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(6.14) is valid for some k. Differentiating (6.14) once more, we obtain 



(i-.^)-^^^M£)=z^+i 



ds^+i 



E^.>71 



j^^x^isy 



J=l 



(1 - z2:gis)r+i 



^,-^^+iV-. (r_+j)(zV^)^-+i 



(1 - z^Jgis^r+i+l 



fc+1 



-'E 






where we have set 



Cj,k+1 '■= < 



Cl,k, J = 1, 

jcj,k + {r + j- l)cj_i,fc, 2 < j <k, 
{r + k)ck,k, j = k + l. 



Hence, by induction, formula (6.14) is valid for all k. 
Now, by (6.14) we have 



mfc(x) = i~ 



dV(s) 



ds^ 



.=0 ^§"^'^l-z-)^-(l-z-')\§'^''' 



since < z^ < 1. Hence, inequalities (6.12) hold with Cfc = c^^ki Ck 

Ek 
j = l Cj,fc. 

6.3. Asymptotics of the moment sums. 
Lemma 6.4. For any k £N and 9 > 0, 



D 



(6.15) 



^.' ' 1 



^dx 



x&X 



(1 - Z^) 



x'lfc 



n 



(fc+2)/3 



n — )■ oo. 



Proof. Using (3.3), by Lemma 3.2 we have 



^ex 



(6.16) 



-e{a,x) 



< 



(j^-{9/2)(a,x) Q^-{e/2){a,x) 



where ao '■= minjai, 02}- On the other hand, 



< 



k \ \k 



^Ox 



(6.17) 



9{a,x) 



^a,x) 



3{a,x) 



> 



{l-z^'f (l-e-<"'^))'^ ~ (a,j;) 



> 
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Since oq ^ \n\~^'^ and \a\ x |n|~^'^, from (6.16) and (6.17) we see that for 
the proof of (6.15) it remains to show 

(6.18) Y. 



g-(a,x) ^ I ,2/3 n^OO. 



Using the Mobius inversion formula (3.10), similarly as in Sections 3 and 
4 we obtain 



i:-- 


-{»,x) ^ Y^ ^(^) Y^ Q-m{a,x} 


xex 


m=l xGZ2_\{0} 




= > Li(m) — — ] 




°° g-moi _|_ g-mQ2 g-m(Qi+Q2) 




■m=l ^ '^ ' 



(6.19) 

By Lemma 3.2, the general term of the series (6.19) is 0{ol{ a^ m~^) (uni- 
formly in m). Hence, by Lebesgue's dominated convergence theorem, the 
right-hand side of (6.19) is asymptotically equivalent to 



^ oo 



fi{m) 






\n 



2/3 



and (6.18) follows. D 

Lemma 6.5. For any A; G N, 

V] |x|'^mfc(x) X |np+^)/^, n -^ oo. 

x£X 

Proof. Readily follows from the estimates (6.12) and Lemma 6.4. D 
Lemma 6.6. For any integer k > 2, 

y^ |x|^/Ufc(x) X |n|(^+^)/^, n-^oo. 

x&X 

Proof. An upper bound follows (for all A; > 1) from the inequality (6.11) 
and Lemma 6.5. On the other hand, by Lyapunov's inequality and formula 
(6.1), for any k > 2 we have 

and a lower bound follows by Lemma 6.4. D 
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Lemma 6.7. For each k £N and j = 1,2 

E:[^,-El{Cj)]''=0{\n\''/'), n^oo. 

Proof. Let j = 1 (the case j = 2 is considered similarly). Using the 
notation (6.9), we obtain, by the multinomial expansion, 

/ \ 2fc 

,2. ^rl^ .A 



2k e 

(6.20) =Y, E ^'^i.--'^^ E 11('^\)''e:[uo{x')'^], 

f-=l ki,...,ke>l, {x\...,xi}cX i=i^ 

ki-\ yki = 2k 

where Cfc^^,,,^fc^ are combinatorial coefficients accounting for the number of 
identical terms in the expansion. Using that E'^[uq{x)] = 0, we can assume 
that ki > 2 for all i = 1, . . . ,^. Since ki + ■ ■ ■ + k^ = 2k, this implies that 
i < k. Hence, recalling the notation (6.10) and using Lemma 6.6, we see 
that the internal sum in (6.20) (over {x^, . . . ,x^} C X) is bounded by 

i i 

E J{\A''^^kA^l<J{Y.\'^\''^^kA^) 

{x\...,xi}CX i=l i=lx&X 



0(1) n hi 



^j(fc,+2)/3 
i=l 

= 0(l)-|n|2(^+^)/3^0(|n|4'=/3), 

and the lemma is proved. D 

7. Local limit theorem. As was explained in the Introduction (see 
Section 1.3), the role of a local limit theorem in our approach is to yield 
the asymptotics of the probability QziC = n} = Q^IIn) appearing in the 
representation of the measure P^ as a conditional distribution, Pn{A) = 
Ql{A\nn) = Ql[A)/Ql{nn), Ac Ha (see (2.16)). 

7.1. Statement of the theorem. As before, we denote a^ := E1[^), K^ := 
Cov(,^,,^) = -EJ(^ — a^)^(C — o^)) where the random vector ^ = (^1,^2) is 
defined in (2.2). From (6.2), it is easy to see (e.g., using the Cauchy-Schwarz 
inequality together with the characterization of the equality case) that the 
matrix K^ is positive definite; in particular, det K^ > and hence K^ is 
invertible. Let Vz = Kz be the (unique) square root of the matrix K~^ 
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(see, e.g., [4, Ch. 6, §5, pp. 93-94]), that is, a symmetric, positive definite 
matrix such that V^ = K~^. 

Denote by /o,/(0 the density of a standard two-dimensional normal dis- 
tribution AA(0, /) (with zero mean and identity covariance matrix), 

/o,/(^) = ;^e-l^l'/2, xGM2. 

ZTT 

Then the density of the normal distribution J\f{az,Kz) (with mean a^ and 
covariance matrix K^) is given by 

(7.1) fa^,KAx) = (.detK,)-^/^foj{{x-a,)V,), x€M.\ 
With these notations, we can now state our local limit theorem. 
Theorem 7.1. Uniformly in m gZ"^, 

(7.2) Q-{e = m} = /,,,^,(m)+0(|7i|-5/3), n^oo. 

Remark 7.1. Theorem 7.1 is a two-dimensional local central limit theo- 
rem for the sum ^ = ^.j.^;^ xi'ix) with independent terms whose distribution 
depends on a large parameter n = (ni, 77-2); however, the summation scheme 
is rather different from the classic one, since the number of non-vanishing 
terms is not fixed in advance and, moreover, the summands actually involved 
in the sum are determined by sampling. 

One implication of Theorem 7.1 will be particularly useful. 

Corollary 7.2. As n ^- 00, 

1/3 

(7.3) Qr{^ = n}^L_^{n,n2r'/\ 



where k = (C(3)/C(2)) 



1/3 



Before proving the theorem, we have to make some (quite lengthy) tech- 
nical preparations, collected below in Sections 7.2-7.4. 

7.2. Lemmas about the matrix norm. The matrix norm induced by the 
Euclidean vector norm | • | is defined by \\A\\ := sup|2.|=]^ \xA\. It is well known 
that for a (real) square matrix A its norm is given by 



(7.4) Pll = VX{A^), 
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where A(-) is the spectral radius of a matrix, defined to be the largest mod- 
ulus of its eigenvalues (see, e.g., [25, §6.3, pp. 210-211]). 

We need some general facts about the matrix norm || • ||. Even though 
they are mostly well-known, specific references are not easy to find (cf., e.g., 
[4, 25, 19]). For the reader's convenience, we give neat proofs of the lemmas 
below based on the spectral characterization (7.4). 

Lemma 7.3 (cf. [17, §22, Theorem 4, p. 40]). If A is a real matrix then 

II /IT /til _ II 4||2 

Proof. The matrix A^A is symmetric and non-negative definite, hence, 
using (7.4), we obtain \\A'^A\\ = X{A^A) = ||^||^, as claimed. D 

Lemma 7.4 (cf. [19, §5.6, Problem 23, hints (2,5) and (5,2), pp. 313-314]). 
If A = {aij) is a real d x d matrix, then 



d d 

2 



(7.5) i E 4 < Mf < E - 



Proof. Note that Yl'i,j=i^'ij = tT{A'^A), where tr(-) denotes the trace, 
and furthermore 

X{A'^A) < tT{A'^A) <d-X{A'^A). 

Since X{A'^A) = \\A\\^ by (7.4), this implies (7.5). D 

The following simple fact pertaining to dimension d = 2 seems to be less 
known. 

Lemma 7.5. Let A be a symmetric 2x2 matrix with detA ^ 0. Then 

\\A\\ 



(7.6) 11^- 



|detA| 



Proof. Let Ai and A2 (IA2I > |Ai| > 0) be the eigenvalues of A, then 
|detj4| = [All • IA2I and, according to (7.4), 



||A|| = yA(l2) = IA2I, 11^-^ = ^A((^"i)2) = |Ai|"\ 
which makes the equality (7.6) obvious. D 
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7.3. Estimates for the covariance matrix. In this section, we collect some 
information about the asymptotic behavior of the matrix K^ = Cov(,^,,^). 
The next lemma is a direct consequence of Theorem 6.1. 

Lemma 7.6. As n ^ oo, 

3(nin2)4/3 



detK, 



1/2 



Let us now estimate the norms of the matrices K^ and Vz = K 
Lemma 7.7. As n — > co, one has \\Kz\\ x |n|^/'^. 
Proof. Lemma 7.4 and Theorem 6.1 imply 

2 

\\Kzf X Y, Kz{i,jf X {nin2f^ x \nf' {n ^ oo), 
«,i=i 

and the required estimate follows. D 

Lemma 7.8. For the matrix Vz = Kz , one has \\Vz\\ x |n|~^''^ as 
n — )• oo. 

Proof. Using Lemmas 7.3 and 7.5 we have 

IIT/ ||2 _ ||T/2|| _ II Ti'-lll _ 11-^^11 

\m -\\Vz\\-\\Kz 11-^^^, 

and an application of Lemmas 7.6 and 7.7 completes the proof. D 

We also need to estimate the so-called Lyapunov coefficient 

(7.7) L, :=||14f ^|x|V3(x), 

where fi^ix) = El\iyo{x)^\ (see (6.10)). 

Lemma 7.9. As n — ^ oo, one has Lz x |n|~^/^. 

Proof. Follows from (7.7) using Lemmas 7.8 and 6.6 (with k = 3). D 
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7.4. Estimates of the characteristic functions. Recall from Section 2.1 
that, with respect to the measure Q^, the random variables {i'{x)}xex are 
independent and have negative binomial distribution with parameters r and 
p = 1 — z^ . In particular, i/{x) has the characteristic function (see (6.13)) 

(7.8) ^.(.)(.) := Elie-'^i^^) = ^^_~fjy , s G M, 

hence the characteristic function ^p^{X) := -E^e*' '^') of the vector ^ = 
Y^x&x^^i^) is given by 

nw = n ^K^)((A,-» = n a^'~/JLY ' ^^^"■ 

Lemma 7.10. Let '^vq(x){^) ^^ ^he characteristic function of the random 
variable 1^0(2^) = i^ix) — El[i'{x)]. Then 

(7.9) |v?^o(x)(s)l <exp{-i/i2(x)sV|//3(a:)|s|3}, seR, 
where fj.k{x) = l^lWo{x)''\ {see (6.10)). 

Proof. Let a random variable ui{x) be independent of 1^0(2;) and have 
the same distribution, and set i>(x) := voi^) — i^i(2;)- Note that El[u{x)] = 
and Var[z>(x)] = 2Var[z^(x)] = 2^2ix). We also have the inequality 

El\i){xf\ < 4:El\uo{xf\ = 4/23 (x) 

(see [5, Lemma 8.8, pp. 66-67]). Hence, by Taylor's formula, the character- 
istic function of i}{x) can be represented in the form 

(7.10) V>u{x){s) = 1 - ti2{x)s^ + |0/U3(2;)S^ 

where \6\ < 1. Now, using the elementary inequality \y\ < e^^ ~i)/2 and the 
fact that 9J^(^.)(s) = \iPyo(x){s)\ , we get 

\^uo{x){s)\ < exp|i(|¥?^(,(^)(s)|^ ~ ^)} " exp|i(v9£)(^)(s) - l)|, 

and the lemma follows by (7.10). D 

The characteristic function of the vector ^0 '■= S, — Oz = Ylxex ^'^oi^) i^ 
given by 

(7.11) ^^,iX) := Er(e'<^'«o)) = n E,^(e^<^'-)'^oW) = J] V'.oW((A,x)). 

xGX x&X 
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Lemma 7.11. If V^ = K^^^'^ then for all A G M^ 

(7.12) \v^oi^V,)\ < exp{-i|A|2 + iL,|A|3}. 
Proof. Using (7.11) and (7.9), we obtain 

(7.13) |^^„(A14)| <exp|-i j;(Ay.,x)V2(x) + ^^KAF.,x)|V3(x)|. 

The first sum in (7.13) is evaluated exactly as 

V(AT4,x)V2(x) = Var(Ay,,0 = XV,CoY{^,0{XV,y 

(7.14) ,^ 

since Cov(,^,^) = K^ = V~'^. For the second sum in (7.13), by the Cauchy- 
Schwarz inequality and on account of (7.7) we have 

(7.15) E KAV^.,^>lV3(x) < \Xf \\V,f Y. l^lV3(x) = |A|%. 
xex xgx 

Now, substituting (7.14), (7.15) into (7.13), we get (7.12). D 

Lemma 7.12. // |A| < L~^ then 

(7.16) |<^g„(Ay.) - e-l^l^/2| < 16L,|A|3e-l^lV6_ 

Proof. Let us first suppose that ^L^^-^^ < |A| < L-^ Then i < L^\Xf < 
|Ap, so (7.12) implies \ip^^{XV,)\ < e-l^l'/^. Hence, 

|995„(AV;) -e"l^l'/2| < |^^,^(Ay,)| +e-l^l'/2 

<2e-l^l'/6<i6L,|A|3e"l^l'/6, 

in accord with (7.16). 

Suppose now that |A| < ^Lz . Taylor's formula implies 

(7.17) 'PuQix){s) - 1 = -^f^2{x)s'^ + le^fX3{x)s^, 
where \dx\ < 1. By Lyapunov's inequality, /i2(a^) < f^six)"^'^, so 

(7.18) \^uo{x){s) - 1| < i|s|V3(a:)'/' + llsfM^)- 
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For s = {XVz,x), we have from (7.15) 

(7.19) KAF.,x)|/i3(x)i/3<Ly3|A|<i, 
and so (7.18) yields 

(7.20) |v5,„(,.)((Ay„x)) - 1| < 1 . 1 + 1 . 1 < 1 . 

Similarly, using the elementary inequality (a + 6)^ < 2(a^ + 6^), from (7.18) 
we obtain 

\^M^)(^) - l|^ < i (|s|/i3(x)^/3 + g|s|V3(a;)) \sffi3{x), 

whence, in view of (7.19) and a general bound |(AT4)2;)| < |A| • \\Vz\\ ■ \x\, it 
follows that 



(7.21) 



k.o(x)((AV^.,x)) - l|2 < 1 (1 + 1 . i)|A|3||F,f |x|V3(^) 



<l\Xf\\Vzf\xffi3{x). 



Consider the function ln(l + y) of complex variable y, choosing the prin- 
cipal branch of the logarithm (i.e., such that In 1 = 0). Taylor's expansion 
implies ln(l + y) = y + Oy^ for \y\ < i, where \e\ < 1. By (7.17), (7.20) and 

(7.21) this yields 

ln(^,„(,)((Ay„x))=-i(Ay„x)V2(x) + i0~.|A|3||y,f|x|V3(x), 

where \9x\ < 1. Substituting this into (7.11), due to (7.14) and (7.15) we 
obtain 

lnv9^„(A14) = ^ln(/.,„((AT4,x)) = - 1 |A|2 + i^iL,|A|=^ (l^il < 1). 

Using the elementary inequality |e^ — 1| < |y|e'^', which holds for any y G C, 
we have 

< e-l^lV2 . iL,|A| V^I"l'/2 < e-l"l'/2L,|A|3, 
and the proof is completed. D 

Lemma 7.13. For all A G M^, 

(7.22) |(^5,(A)|<exp{-irJ„(A)}, 
where 

(7.23) J„(A) := Y^ e-<"'^>(l - cos(A,x)). 
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Proof. According to (7.11), we have 

(7.24) |<^5,(A)| = |(^g(A)| = expJ ^ ln|(p,(,)((A, x))\ I. 

Using (7.8), for any s G M we can write 

r ll — z^\ T ( ll — z^\ \ 

ln|^.(.)(.)| = - In ^3^-^ < - y^^^—^. - Ij 

r2:^(l — coss) rz^(l — coss) 



Utilizing this estimate under the sum in (7.24) (with s = (A, x)) and recall- 
ing the notation (3.3), we arrive at (7.22). D 

7.5. Proof of Theorem 7.1 and Corollary 7.2. Let us first deduce the 
corollary from the theorem. 

Proof of Corollary 7.2. According to Theorem 5.1, a^ := El{^) = 
n + o(|np^). Together with Lemma 7.8 this implies 

|(n-a^)V;| < In-Ozl • ||14|| 

= o(|n|2/3)0(|n|-2/3) =o(l). 

Hence, by Lemma 7.6 we get 

/a.,i..(n) = :^(deti^.)-V2e-l(--)V'^lV2 



''' (nin2)-2/=^(l + o(l)). 



r ' K 



2\/3 7r 
and (7.3) follows from (7.2). D 

Proof of Theorem 7.1. By definition, the characteristic function of 
the random vector ^o = S, ~ o-z is given by the Fourier series 

^5o(A) := E,^(e*<^'«o)) = ^ Ql{^ = m} e^^^'™--^^), A G R^ 

hence the Fourier coefficients are expressed as 

(7.25) Ql{^ = m} = ^l^ e-*<^' — '^^^ (^^^ ( A) dA, m G Z^ , 
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where T^ := {A = (Ai,A2) G M^ : |Ai| < vr, IA2I < tt}. On the other hand, 
the characteristic function corresponding to the normal probabihty density 
fa,,KAx) (see (7.1)) is given by 

^a.,i..(A) = e^<^'->-l^^^''l'/2, AGM2, 

so by the Fourier inversion formula 

(7.26) fa.,KM) = A / e-^^'^-'^^^-l^^'l'/^dA, m G Z^. 

Note that if lAy^""*^! < Lj^ then, according to Lemmas 7.8 and 7.9, 
|A| < \XV-'\ ■ \\V,\\ < L-^\\V,\\ = 0(|nr^/3) ^ ^^^^^ 

which implies that A G T^. Using this observation and subtracting (7.26) 
from (7.25), we get, uniformly in m G Z^, 

(7.27) \Ql{^ = m}- fa.,KAm)\ <h+h + h, 
where 

h:=^J ^ J,,,,(A)-e-l^^^-^l^/^|dA, 

h:=^J e-l^^^~^l^/^dA, 



h--=^,l l^eo(A)|dA. 



47r2 



T2n{A:|Ay,->LJ^} 



By the substitution \ = yVz, the integral /i is reduced to 

IdetT^I [ I , ^,. u,i2/ 



47r2 /|„|<r- 

(7.28) •^lyl<^. 



/ |(^^^(yy,)-e"IJ'l'/2|dy 

J|?/|<L7^ 



0(l)(deti^,)-i/2L, /" |y|Vl^l'/6dy = 0(|nr5/3)^ 

on account of Lemmas 7.6, 7.9, and 7.12. Similarly, again putting \ = yVz 
and passing to the polar coordinates, we get, due to Lemmas 7.6 and 7.9, 

Idety.l /■- _|^|2/2 ,, , 
h = — 7, / \y\ e '^' ' a\y\ 



(7.29) 27r 7i-i 

= 0(|n|-^/3)^-L-/2^^(|^|-5/3)_ 
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Estimation of I3 is the main part of the proof. Using Lemma 7.13, we 
obtain 



(7.30) h = 0{l) / e--^"WdA, 

JT^n{\xvr'^\>L-^} 

where Ja(A) is given by (7.23). The condition jAF^"^] > L^^ implies that 
|A| > rj\a\ for a suitable (smaU enough) constant rj > and hence 

max{|Ai|/ai, |A2|/a2} > 7], 

for otherwise from Lemmas 7.7, 7.9 and 7.12 it would follow 

1 < L^\XV-'^\ < L,r]\a\ ■ \\K^\\^/^ = 0{r]) -^ as r/ i 0. 

Hence, the estimate (7.30) is reduced to 

(7.31) I3 = 0(1) ( / + f I e-^'^W dA. 

To estimate the first integral in (7.31), by keeping in the sum (7.23) only 
pairs of the form x = {xi, 1), xi £ Z+, we obtain 

00 
e"V„(A) > Y^ e-"i^i (1 - Ke*(^i^i+^2)) 



a;i=0 
1 


-K 


r Qi^2 


1 - e-^^i 

1 


1 



(7-32) > ^ ,^ ^., , , 

because Ku < \u\ for any u E C Since -qai < |Ai| < vr, we have 

|^_g-ai+iAij > j-^_g-ai+ir,aij _ ^^(^^ ^ ^2)1/2 (^^ ^q)^ 

Substituting this estimate into (7.32), we conclude that Ja(A) is asymptoti- 
cally bounded from below by C(r/)a^ X |n|^''^, uniformly in ryai < |Ai| < vr. 
Thus, the first integral in (7.31) is bounded by 

0(1) exp (-const • \n\^/^) = o{\n\^^l'^). 

Similarly, the second integral in (7.31) (where IA2I > 7702) is estimated by 
reducing summation in (7.23) to that over x = (1,2:2) only. As a result, we 
obtain that I^, = o(|n|~^/^). Substituting this estimate together with (7.28) 
and (7.29) into (7.27), we get (7.2), and so the theorem is proved. D 
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8. Proof of the limit shape results. Recall the notation (see (4.1), 
(4.2)) ^(t) = J2xex{t)^'^i^)^ where X{t) = {x £ X : xa/xi < t(n2/ni)}, 
t G [0,oo]. As stated at the beginning of Section 4, the tangential parame- 
terization of the scaled polygonal line Fn = Sn (r) is given by 

(8.1) Ut) := Snim) = K'6(i), n^^Ut)), t G [0,oo], 

whereas the limit shape 7* determined by equation (1.3) is parameterized 
by the vector- function g*{t) = {gl{t),g2{t)) defined in (4.4) (see more details 
in the Appendix, Section A.l). 

The goal of this section is to use the preparatory results obtained so far 
and prove the uniform convergence of random paths Cn(') to the limit g*{-) 
in probability with respect to both Q^ (Section 8.1) and P^ (Section 8.2). 
Let us point out that, in view of (8.1), Theorems 8.1 and 8.2 below can 
be easily reformulated (cf. Theorem 1.1 stated in the Introduction) using 
the tangential distance d'j-{rn,j*) = supo<t<oo|Cn(0 ~9*{'t)\ (see (1.9); cf. 
general definition (A. 3) in Section A.l below). 

8.1. Limit shape under Q!]. Let us first establish the universality of the 
limit shape under the measures Q^, which in conjunction with the next The- 
orem 8.2 illustrates the asymptotic "equivalence" of the probability spaces 
(il, Q^) and (77„,P„^). 

Theorem 8.1. For each e > 0, 

hm Ql\ sup \n-'Cj{t) - g*{t)\ < e] = 1 {j = 1,2). 

Proof. By Theorems 4.1 and 4.2, the expectation of the random process 
n~ S,j{t) uniformly converges to g1{t) as ?i — )• 00. Therefore, we only need 
to check that for each e > 

hm Qll sup n-%{t) - El[C,{t)]\ > e] = 0. 



Note that the random process Coj(i) '■= Ci(^) ~ ^liCjit)] has independent 
increments and zero mean, hence it is a martingale with respect to the 
natural filtration J^t '■= <7{z^(x), x G X{t)}, t G [0, 00]. From the definition of 
S,j{t) (see (4.2)), it is also clear that ^oj{t) is a cadlag process, that is, its paths 
are everywhere right-continuous and have left limits. Therefore, applying the 
Kolmogorov-Doob submartingale inequality (see, e.g., [39, Corollary 2.1, 
p. 14]) and using Theorem 6.1, we obtain 

Q:\ sup |eo,(t)|>enA<^^^ = 0(|n|-2/3)^0, 

L0<t<oo J {^^j) 
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and the theorem is proved. D 

8.2. Limit shape under P^'. We are finahy ready to prove our main re- 
sult about the universality of the limit shape under the measures P^ (cf. 
Theorem 1.1). 

Theorem 8.2. For any e > 0, 

hm P;;\ sup \nj'C,it) - g*it)\ < e] = 1 (j = 1, 2). 

Proof. Similarly as in the proof of Theorem 8.1, the claim of the theorem 
is reduced to the limit 



lim Pl^i sup koi(i)| > C"-? r = 0' 

n^oo [o<t<oo J 



where Coj{t) = Cj{t) - El[Cj{t)]. Using (2.9) we get 

'3z{supo<t<ooKoj(i)| > ^^j} 



i.2) P:{ sup \(oj{t)\ >enj]< 

L0<t<oo J 



QliC = n} 



Applying the Kolmogorov-Doob submartingale inequality and using Lemma 
6.7 (with A: = 3), we obtain 

L0<t<oo J [^^j) 

On the other hand, by Corollary 7.2 

Q;{^ = n}x(nin2)-2/3x|n|-4/3. 

In view of these estimates, the right-hand side of (8.2) is dominated by a 
quantity of order of 0(|n|~^/^) — )• 0, and the theorem is proved. D 

APPENDIX 

A.l. Tangential distance between convex paths. Let Qq be the 

space of paths in M^ starting from the origin and such that each path 7 E ^g 
is continuous, piecewise C^-smooth (i.e., everywhere except a finite set), 
bounded and convex, and, furthermore, its tangent slope (where it exists) 
is non- negative, including the possible value -|-oo. Convexity implies that 
the slope is non-decreasing as a function of the natural parameter (i.e., the 
length along the path measured from the origin). 
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For 7 G Qq, let ^^(i) = {gi{t),g2{t)) denote the right endpoint of the 
(closure of the) part of 7 where the tangent slope does not exceed t G 
[0,00]. Note that the functions xi = gi{t), X2 = g2{t) are cadlag (i.e., right- 
continuous with left limits), and 

(A.l) ^ = 474 = ^- 

More precisely, equation (A.l) holds at points where the tangent exists and 
its slope is strictly growing; corners on 7 correspond to intervals where both 
functions gi and 52 are constant, whereas flat (straight line) pieces on 7 lead 
to simultaneous jumps of gi and 52- 

The canonical limit shape curve 7* (see (1.3)) is determined by the para- 
metric equations (cf. (1.8), (4.4)) 

^2 _|_ 2t t^ 

(A.2) XI = gl{t) = ——-^ , X2 = gl{t) = ——-^, tG[0,oo]. 

Indeed, it can be readily seen that the functions (A.2) satisfy the Cartesian 
equation (1.3) for 7*; moreover, it is easy to check that t in equations (A.2) 
is the tangential parameter: 



dglit) 


2 dgl 


2t 


dt 

and hence (cf. (A.l)) 


{l + tf dt 

dx2 _ dgl/dt _ 
dxi dgl/dt 


(l + t)3' 



The tangential distance dj- between paths in Qq is defined as follows: 

(A.3) (i7-(7i,72) := sup |57i(*) - fi'72(*)li 7i,72 G ^o- 

0<i<oo 

Lemma A.l. The Hausdorff distance d-^ defined in (1.2) is dominated 
by the tangential distance d-j-: 

(A.4) c?h(7i,72) < c^r(7i:72), 7i,72e^o- 

Proof. First of all, note that any path j £ Qo can be approximated, 
simultaneously in metrics d-^ and df, by polygonal lines 7"^ (for instance, 
by inscribing polygonal lines with refined edges in the arc 7) so that 

lim dw(7,7") = 0, lim dr(7,7"') = 0. 
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This reduces inequality (A. 4) to the case where 71, 72 are polygonal lines. 
Moreover, by symmetry it suffices to show that 

(A. 5) maxmin |x — y| < d7-(7i,72). 

a;e7i J/G72 

Note that if a point x € 71 can be represented as x = g^-^ (to) with some 
to (i-e., X is a vertex of 71), then 

min|3;-y| = min|5^^(to) - y\ < Ig^A^o) - d-yii^o)] < t^r(7i,72), 
yG72 y&-y2 

and inequality (A. 5) follows. 

Suppose now that x G 71 lies on an edge — say of slope t* — between two 
consecutive vertices g^-^{t^) and g^^t*), then x = Og^^t^) + (1 — 9)g^At*) 
with some G (0, 1) and 

min|x -y\ = m.ui\e g^-^{t^) + (1 - 9)g.yAt*) - y\ 
4/672 3/G72 

< \eg^,{Q + {1 - 6)g^,{f) - g^,{t*)\ 

< dlg^AQ - 9'f2{t*)\ + (1 - G)\gr,{t*) - gj2{i*)\ 
(A.6) < e\g^,{U) - g^,{t*)\ + (1 - 9) dr{-f 1,72). 

Note that for all t G [t^,t*) we have g^/^t) = g-y^tt,), hence 

(A.7) I571 (**)- 572 Wl = l57i (0-572 Wl <rfr(7i, 72) {U<t<t*). 

If 572 (t) is continuous at t = t*, then inequality (A.7) extends to t = t*: 

\g-yi{t*) - gj2{i*)\ < rfr(7i,72)- 

Substituting this inequality into the right-hand side of (A.6), we see that 
miuyg^j \x — y\ < df{'yi,^2), which implies (A. 5). 

If 5^2 (t* — 0) 7^ 572 (^*) then t* coincides with the slope of some edge on 72 
(with the endpoints, say, g^^iK) ^iid g')2{^*))-> which is thus parallel to the 
edge on 71 where the point x lies (i.e., with the endpoints 571 (t*), 57^ (t*)). 
Setting s* := max{ t'^,t*} < t*, we have g^^t^) = g^^s^), 5'72(it) = 5'72(s*)- 

To complete the proof, it remains to observe that the shortest distance 
from a point on a base of a trapezoid to the opposite base does not exceed 
the maximum length of the two lateral sides. Hence, 

min|x-y| < minjlx-yl : y G [572('5*),5'72(^*)]} 
yG72 

< max{|5f7^(s*) - 5'72('S*)M57i(**) -572(**)l} 

< (ir(7i,72), 

and the bound (A. 5) follows. D 
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Remark A.l. Note, however, that the metrics d-^ and dj- are not equiva- 
lent. For instance, if 7 G ^o is a smooth, strictly convex curve with curvature 
bounded below by a constant xq > 0, then for an inscribed polygonal line 
r^ with edges of length no more than e > 0, its tangential distance from 7 
is of the order of e, but the Hausdorff distance is of the order of e^: 

dr(n,7)xe, dK(r„7)xxoe2 (e ^ 0). 

Moreover, in the degenerate case where the curvature may vanish, the differ- 
ence between the two metrics may be even more dramatic. For instance, it is 
possible that two polygonal lines are close to each other in the Hausdorff dis- 
tance while their tangential distance is quite large. For an example, consider 
two straight line segments A 5-^2 C M^ of the same (large) length L, both 
starting from the origin and with very close slopes, so that the Euclidean dis- 
tance 5 between their right endpoints is small; then (i-^(/^i, 12) < <^ whereas 
dr(A,r2) = L. 

A. 2. Total variation distance between P^ and P^. Note that if 
probability measures Pn and Pn on the polygonal space Tin are asymptoti- 
cally close to each other in total variation {TV), that is, \\Pn — PuWtv -^ 
as n — >• 00, where 

\\Pn-Pn\\TV-= sup |P„(^) - P„(A)| , 

then the problem of universality of the limit shape is resolved in a trivial way, 
in that if a limit shape 7* exists under P„ then the same curve 7* provides the 
limit shape under P„. Indeed, assuming that the event Afr = {d{Tn,j*) > e} 
satisfies PniAgr) — )• as n — )• (X), we have 



^e; 



PniAs) < Pn{Ae) + |P„(^e) " ^n(Ae , 

<Pn{Ae) + \\Pn-Pn\\TV^^ (n^Oo). 

However, the family {Pn}, defined by formula (2.10) with the coefficients 
(2.11), is not close to Pn in total variation, at least uniformly in r E (0, 00). 

Theorem A. 2. For every fixed n, the limiting distance in total variation 
between Pn and P^, as r — t- or r —t- 00, is given by 

(A.8) lim \\P: - P^Wtv = I 
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Proof. To obtain a lower bound for ||P^ — P^H^y in the case r — )• oo, 
consider the polygonal line F* G 77„ consisting of two edges, horizontal 
(from the origin to (ni,0)) and vertical (from (ni,0) to n = (^1,722)). The 
corresponding configuration vp* is determined by the conditions z//^*(l,0) = 
^ii z^r*(0, 1) = n2 and vr*{x) = otherwise. Note that Wj, ~ r /k\ as 
r — )• 00 (see (2.11)), hence U'{r) = 0{r^'^^, where Np := '^x(^x^r{x) is 
the total number of integer points on F (excluding the origin). We have 
Nr* = ni + ns, so V'iF*) = V^^U^^ ~ r^i+^VK'^2!) as r ^ 00. On the 
other hand, for any F e n^ {F ^ F*) one has 6^(r) = o(r"i+"2). Indeed, 
for X S ^ we have xi + X2 > 1 and, moreover, xi + a;2 > 1 unless x = (1,0) 
or X = (0, 1). Hence, for any T G 77„ (r / F*), 

Nr = ^ J^rix) < ^ (xi + X2) i^rix) = ni + n2, 

so that Nr < ni+n2 and b^'{F) = 0{r^'^) = o(r"i+"2) as r — )• 00. Therefore, 
from (2.10) we get 

ff(r*\ „ni+n2 

n^ ^ &'"(^*) + Er^r-^''(-^) r"i+"2+o(r"i+"2) ^ ^' 

and it follows that 

1 



(A.9) i|p; - p^]\tv > Kin - Pnini ^ 1 



#(^n) 



00 



For the case r — )• 0, consider the polygonal line F^ G 77„ consisting of one 
edge leading from the origin to n = (ni,n2). That is, z^r*('^/^n) = ^n and 
z^(x) = otherwise, where /c„ := gcd(ni,n2). Clearly, f'^F^,) = b'^. = r/kn, 
while for any other F £ Fin (i.e., with more than one edge), by (2.11) we 
have b^{F) = 0{r'^) as r — )• 0. Therefore, according to (2.10), 

P'(r ) = ^'^^^*^ = - yl (r ^0) 

and similarly to (A.9) we obtain 

(A.IO) \\P: - P^Wtv > Kin) - P'AF.)\ ^ 1 - -77^ (r ^ 0). 

The upper bound for \\P^ — P^\\tv (uniform in r) follows from the known 
fact (see [10, p. 472]) that the total variation distance can be expressed in 
terms of a certain Vasershtein (-Kantorovich-Rubinstein, cf. [36]) distance: 

ll^n-^nllTy=infE[£,(X,y)]. 

Ji. , Y 
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Here the infimum is taken over all pairs of random elements X and Y defined 
on a common probability space (fl, J-, P) with values in i7„, and the marginal 
distributions P^ and P^, respectively; E denotes expectation with respect 
to the probability measure P, and the function q{-, •) on iT„ x -^n is such 
that g{r, r') = lifr^ r' and Q{r, r') = if r = r' (therefore defining a 
discrete metric in Tin). Choosing X and Y so that they are independent of 
each other, we obtain 

\\P: - PuWtv < ^[QiX,Y)] = 1-P{X = Y} 

= 1- ^ P{x = r,Y = r} = i- Y, Kin-P'nin 

Combining this estimate with (A. 9) and (A. 10), we obtain (A. 8). D 

In the limit n — )■ oo. Theorem A. 2 yields 

lim lim IIP^ - P„^||Ty = 1- 

n— >-oor— >0,oo 

That is to say, in the successive limit r — >• (oo), n — ;• oo, the measures P^ 
and P}^ become singular with respect to each other. 

Moreover, one can show that the distance ||-P^ — P^||tv is not small even 
for a fixed r / 1. To this end, it suffices to find a function on 77^ possessing 
a limiting distribution (possibly degenerate) under each P^, with the limit 
depending on the parameter r. Recalling Remark 2.1, it is natural to seek 
such a function in the form referring to integer points on P £ Pin- Indeed, 
for the statistic Np = '^Zx&x ^r{x) introduced in the proof of Theorem A. 2, 
the following law of large numbers holds (see Bogachev and Zarbaliev [6, 
Theorem 3] and Zarbaliev [40, §1.10]). 

Lemma A. 3. Under Assumption 3.1, for each r G (0, oo) and any e > 0, 



(A.ll) lim Pl{Al) = 1, where A^ :-- 



< e 



{nin2y^^ /^^ 

The result (A.ll) implies that, for any r ^ 1 and all e > small enough, 
\\P: - P^Wtv > KiAl) - P^{Al)\ ^ |1 - 0| = 1 (n ^ oo), 
and we arrive at the following result. 
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Theorem A. 4. For every fixed r ^ \, we have 

lim \\PI-P},\\tv = 1, 

n— >oo 

hence the measures P'^ and P^ on 11^ are asymptotically singular with re- 
spect to each other as n — t- oo. 
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